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Abstract 

We propose an alternative implementation of the curvaton mechanism for generating the 
curvature perturbations which does not rely on a late decaying scalar decoupled from 
inflation dynamics. In our mechanism the supersymmetric Higgs scalars are coupled to 
the inflaton in a hybrid inflation model, and this allows the conversion of the isocurvature 
perturbations of the Higgs fields to the observed curvature perturbations responsible for 
large scale structure to take place during reheating. We discuss an explicit model which 
realises this mechanism in which the p term in the Higgs superpotential is generated 
after inflation by the vacuum expectation value of a singlet field. The main prediction 
of the model is that the spectral index should deviate significantly from unity, \n — 1| ~ 
0.1. We also expect relic isocurvature perturbations in neutralinos and baryons, but no 
significant departures from gaussianity and no observable effects of gravity waves in the 
CMB spectrum. 
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1 Introduction 


According to the inflationary paradigm [1], the presently observed large scale universe originated 
from a very small patch of space which underwent a period of quasi-exponentially accelerated 
expansion known as inflation. In such an inflationary approach, the very largest scales, which 
are now entering the horizon, would have been in causal contact at very early times. According 
to inflation this causal contact would have ceased some 50-60 e-folds before the end of inflation, 
corresponding to so-called horizon exit to distinguish it from the horizon entry as observed in the 
present epoch. Inflation clearly provides an attractive explanation for why the cosmic microwave 
background (CMB) radiation should have the same temperature from points in the sky which 
would have been out of causal contact at the time of last scattering (the horizon problem). It 
also accounts for the observed flatness of the universe hi = 1 (the flatness problem), consistent 
with the CMB data [2-7]. 

Another commonly stated success of inflation [1] is the fact that the observed primordial 
density perturbations, which were first observed by COBE [8] on cosmological scales just entering 
the horizon, and which are supposed to be the seeds of large scale structure, could have originated 
from the quantum fluctuations of the inflaton held, the scalar held which is supposed to be 
responsible for driving inflation. In this scenario the quantum fluctuations of the inflaton field 
during the period of inflation become classical perturbations at horizon exit, giving a primordial 
curvature perturbation which remains constant until the approach of horizon entry [9]. The 
advantage of this scenario is that the prediction for the nearly scale-invariant spectrum depends 
only on the form of the inflaton potential, and is independent of what goes on between the 
end of inflation and horizon entry. The disadvantage is that it provides a strong restriction 
on models of inflation. The price of such simplicity, with one held being responsible for both 
inflation and the primordial curvature perturbation often translates into a severe restriction on 
the parameters of the inflaton potential. This often requires very small values for the couplings 
and/or the masses which apparently renders many such theories unnatural. 

Recently it has been pointed out that in general it is unnecessary for the inflaton field to 
be responsible for generating the curvature perturbation [10,11]. It is possible that the inflaton 
only generates a very small curvature perturbation during the period of inflation, which instead 
may result from the isocurvature perturbations of a curvaton held which subsequently become 
converted into curvature perturbations in the period after inflation, but before horizon entry 
[12,13,10,11]. Isocurvature perturbations simply mean perturbations which do not perturb the 
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total curvature, usually because the curvaton field contributes a very small energy density 
during infiation. In the scenarios presented so far [10,11,14-17], the curvaton is assumed to 
be completely decoupled from inflationary dynamics, and is assumed to be some late-decaying 
scalar which decays before the time that neutrinos become decoupled. 

The reason why the curvaton is assumed to be late-decaying can be understood from the 
following argument [16]. After reheating the total curvature perturbation can be written as, 

77 = (1 - f)n, + fn^ ( 1 ) 


where 


/ = 


3po 


4pr + 3Pa 

and (on unperturbed hypersurfaces on snper-horizon sized scales) 


77, -H 


' 5 Pi 
. Pi 


'5^ 


( 2 ) 


( 3 ) 


where the curvaton density pa and radiation density pr, arising from the decay of the infiaton, 
each satisfy their own energy conservation eqnations and each 77j remains constant on super¬ 
horizon scales. The time evolution of 77 on these scales is then given by its time derivative. 


77 ^-77/(1-/)^ 
where S^r is called the entropy perturbation defined by 


( 4 ) 


Sar ^ -3(77,^ - 77r 


( 5 ) 


The curvaton generates an isocurvatnre pertnrbation because initially pr 3> Po-, and hence 
/ <^ 1, so that from Eq. (1) the curvatnre perturbation is dominated by 77r. However as the 
nniverse expands and the scale factor a increases while the Hubble constant 77 decreases, the 
curvaton with mass m, whose oscillations have effectively been frozen in by the large Hubble 
constant, begins to oscillate and act as matter. After this happens pr decreases as while 
the energy density in the curvaton field p^ has a slower fall-off as a~^. Eventually the curvaton 
energy density pa becomes comparable to the radiation density from the infiaton decay pr, and 
when this happens from Eq. (2) we see that / ~ 1 and from Eqs. (4), (5) this leads to the 
growth of the total curvature perturbation 77 from the isocurvature perturbation 77o- > 77r. This 
mechanism, which allows the curvaton isocurvature perturbations to become converted into the 
total curvatnre perturbation, requires the curvaton scalar to be late-decaying. 
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The main motivation behind the present paper is to propose an alternative curvaton mech¬ 
anism that removes the necessity for having a late-decaying scalar. To this end we make the 
following two observations: 

(1) The first observation is that hybrid inflation [18,19] involves more than one scalar field 
and ends by a phase transition in which the fields involved move from the false, vacuum energy 
dominated potential minimum, towards the global minimum, and start oscillating around it. 
In the process, the vacuum energy gets redistributed among the fields such that their energy 
densities are comparable. Thus any isocurvature perturbation in one of the hybrid inflation 
fields may be converted into curvature perturbations during the on-set of reheating. Note that 
the conversion does not take place until one of the fields decays. Given that the inflaton field is 
a singlet field, we can choose the curvaton field to be a flat direction during inflation made of a 
pair of charged fields, where the gauge interactions imply that the curvaton decays first. 

(2) The second observation is that the supersymmetric standard model provides a natural 
candidate for such a pair of charged scalar fields: the two Higgs doublets. In order to allow a 
flat direction during inflation, we also require that the p term which usually couples the two 
Higgs doublets in the minimal supersymmetric standard model must be set to zero. However, 
providing this obstacle can be overcome, it is natural to explore the attractive possibility that 
such Higgs fields are the seed of large-scale structure in the Universe. The idea would be that 
the inflaton gives a very small curvature perturbation, with the Higgs giving an isocurvature 
perturbation. The coupling of the Higgs scalars to the inflaton then allows the energy densities 
of the inflaton and the Higgs to become of similar magnitude at the on-set of reheating, allowing 
the conversion of isocurvature density perturbations to curvature perturbations. 

In the remainder of paper we shall discuss a supersymmetric hybrid inflation model based 
on these observations. In this model the superpotential includes, in addition to the two Higgs 
doublet superfields Hu and Hd of the supersymmetric standard model, also two gauge singlet 
superfields, one of them 0 playing the role of the inflaton and the other N coupling to the 
two Higgs doublets, effectively generating an effective p term from the coupling NHuHd as 
in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [20]. This happens after 
inflation when it gets a vacuum expectation value (vev) Nq =< N >. The N superfield also 
couples to the inflaton 0 as 0iV^, and therefore acts as a messenger allowing the Higgs fields 
to couple to the inflaton. In previous discussions of this model [20] we assumed that during 
inflation the Higgs and N fields are held at the origin, but here we shall show that an alternative 
inflationary trajectory is possible in which the inflaton 0 as well as the Higgs doublets are slowly 
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rolling during the inflationary epoch. The isocurvature perturbations of the Higgs fields during 
inflation are converted into curvature perturbations during the initial stages of the reheating 
process, when all these fields begin to oscillate with comparable energy densities. 

The layout of the rest of the paper is as follows. The model is presented in section 2, 
where the evolution for the background fields is studied. In section 3 we discuss the nature 
and evolution of the perturbations during the epoch of inflation, and show that it is natural 
for the inflaton to give the dominant curvature perturbation, but still too small to account for 
the COBE value, while the Higgs fields strongly contribute to the isocurvature perturbations. 
In section 4 we describe the evolution of the perturbations during reheating, after the Higgs 
has decayed and the Universe is made of a mixture of matter (the oscillating singlet fields) and 
radiation (Higgs decay products), and show that the isocurvature perturbations are converted 
into curvature perturbations. In section 5 we comment on some of the subtleties involved in 
the transition from inflation to reheating, and the issue of preheating. The predictions of the 
approach are discussed in section 6, and in section 7 we provide a summary. 

2 The model: a new inflationary trajectory 

In this section we revisit the supersymmetric hybrid inflation model based on the superpotential 

[ 20 ]: 

W = \NHuHd - K(t)N^ , (6) 

where N and 0 are singlet superflelds, and are the Higgs superflelds, and A, k are dimen¬ 
sionless couplings. Other cubic terms in the superpotential are forbidding by imposing a global 
U{l)pQ Peccei-Quinn symmetry. The superpotential in Eq. (6) includes a linear superpotential 
for the inflaton field, 0, typical of hybrid inflation, as well as the singlet N coupling to Higgs 
doublets as in the NMSSM. In the original version of this model [20] we assumed that during 
inflation N and Hu, H^ were set to zero. Here we discuss an alternative inflationary trajectory in 
which these fields may take small values away from the origin, consistent with slow roll inflation. 

In order to satisfy the D-flatness condition^, we assume the values of the Higgs doublets 
during inflation to be equal, = H^ = h. Inflation takes place below the susy breaking scale. 
Therefore, including the soft susy breaking masses, rrijsi and rrih-, and trilinears Ax, the 

^Higgs fields are charged under SU{2)l x U{\)y interactions, such that there is what is called a D-term 
contribution in the potential of the form {g^ + g\){Hu — /8. A “D-flat” direction is made of a combination 

of the fields such that the D-term vanishes [21]. 
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potential for the real part of the fields is: 

^ = l/(O) + ^iV^ + «2(0-0+)(0-OiV2 + lmj02 

+\Nh^{-^ — K(j)) + —h^ + + -m\h^ . (7) 

In our previous work we set the Higgs terms to zero, h = 0. Allowing the Higgs terms in Eq. 
(7) we see that the term proportional to XNh"^ can induce a non-zero value for h and N during 
inflation providing that the coefficient of this term is negative ^ — k4> < 0. However, providing 
that the values of the helds N, and are sufficiently small during inflation, as discussed 
later, then inflation is controlled by the inflaton 0, as in the original model [20], and will end 
when 0 reaches one of its critical values 0^ given by: 


K0^ 






( 8 ) 


These critical values correspond to the held dependent mass squared /t^(0 — 00)(0 — (f>~)N‘^ 
changing sign and becoming negative, signalling the end of inhation due to the destabilisation 
of the inhationary trajectory, after which the helds N, 0, Hu and then approach their global 
minimum and acquire their physical vevs Nq, 0o, Vu and Vd- For order of magnitude estimations, 
here on we will take 0^ ~ 00 — 07- 

The parameters of the potential in Eq. (7) are selected by the following physical requirements. 
As discussed in [20] the global PQ symmetry imposed on the superpotential Eq. (6) solves the 
strong CP problem. When the PQ symmetry is broken by the vevs of the singlets, it leads to a 
very light axion in the usual way [22]. The axion scale fa is set by the vevs of the singlets, and 
it is constrained by astrophysical and cosmological observations to be roughly in the window 
10^^ GeV < fa < 10^^ GeV [23,24]. Given that /« ~ 0c ~ 0o ~ -^o ~ 10^^ GeV and the soft 
breaking term are expected to be of the order of 1 TeV, Eq. (8), this leads to a coupling 
constant of the order k ~ 10“^°. The same applies to A ~ k, with fi = XNq ~ ITeH. The 
smallness of A and k will be explained in a companion paper [25]. Therefore, demanding a zero 
cosmological constant at the global minimum requires the height of the potential during inflation 
to be of the order of 

H(O)i/"~^0e~(lOSG'eH), (9) 

with a Hubble parameter of the order of 


H 


—M~ lOMeH, 

v^mp 


( 10 ) 
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where mp = Mp /= 2.4 x 10^® Gel^ is the reduced Planck mass. 

An important condition for inflation is that the inflaton mass (and also rrih) needs to 
be small enough in order to ensure the slow-roll of the inflaton, as determined by the slow roll 
parameters defined in the standard way as [ 1 ]: 


\v\ 



mlV” 

V 


< 1 


( 11 ) 

( 12 ) 


where V'{y'') are the first (second) derivatives of the potential. In this case we hnd the slow 
roll parameters for the 0 field: 

<1, (13) 


m. 


= 




1 2 0 ^ 


(14) 


where it is understood that these expressions are evaluated at some number of e-folds before the 
end of inflation. Eq. (14) is always fulhlled once 77 ^ < 1. This would require an inflaton mass 
of the order of some MeVs at most, say ~ 0.1 A similar constraint will also apply to 
the Higgs soft masses. This is also compatible with the observational constraint on the spectral 
index, n = 0.93 ±0.13 [2], which in terms of the slow-rolling parameters [1], | n — 1 |= 277 ^ — Oe^, 
gives 770 < 0.03. The smallness of the inflaton and Higgs soft masses which are orders of 
magnitude smaller than the typical soft breaking values assumed for the trilinear couplings, is 
explained in a companion paper [25]. 

We have also included a constant vacuum energy contribution V (0) in the potential in Eq. (7) 
whose origin we do not specify. It can be explained once the model is embedded in a supergravity 
(sugra) model [26], where the question of susy breaking can be addressed. It should be pointed 
out however that sugra corrections generically give rise to scalar masses contributions of the order 
of the Hubble parameter during inflation. This is the so-called 77 -problem [27]. Nevertheless, the 
problem can be avoided [19,26,28,29] by a suitable choice of the sugra model, i. e., the Kahler 
potential and/or the superpotential for the scalar inflaton. Our approach is that it is natural 
to have scalar masses of the order during inflation, but rather slightly smaller in order to 
satisfy Eq. (13). This would only require at most a mild tuning of the parameters in the sugra 
Kahler potential. 

We now discuss an inflationary trajectory that will enable the Higgs helds to be non-zero and 
slowly rolling during the inflationary epoch, and hence to acquire an isocurvature perturbation. 
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The requirement that the h, (p fields be slow-rolling during infiation, means that their effective 
field dependent masses must be smaller than 


d^V 

'W 

d^V 

d(j)^ 


~ A2(3h2 + iv2) + 2AiV(4^ 

V2 

~ 


K(j)) -I- m\ < 


(15) 

(16) 


This requirement implies that both the fields N and h must take small values during infiation, 
with 


N h 

^<10-10, ^<10-^ (17) 

0c 0c 


and also, as already remarked, must have soft masses of order an MeV or less. 
The N field-dependent mass is much larger and positive. 


d^V 


~ 2^(4, - 4t){4 - C) ~ 0(TeV2), 


(18) 


where we dropped the small term using Eq. (17). Thus the N field will oscillate with an 

amplitude damped by the expansion, following the evolution equation: 


N + 3HN + 


dV 


0 , 


(19) 


where 


^ (2^2(0-0+)(0-0^) + A2/i2)iV + Ah2(^-K0) 

~ ^^(0)^ +Ah^(^ - K0), (20) 

and we have again dropped the small term in the first line using Eq. (17). The frequency 
of the oscillations is then of the order ~ Kcpc- 

Apparently, as far as N is concerned, we have violated the slow roll conditions. Since we have 
seen that N does not disturb the slow roll of h and (p we should not be concerned about this. 
Nevertheless, we shall see that in some sense N can be regarded as slowly rolling, according to 
the following argument. The typical time scale during infiation is set by the rate of expansion 
H, but the oscillations of N are much faster than that. On the Hubble time scale as far as 
the motion of N is concerned all that we can see is the average effect of the oscillations. The 
motion of N is then given by the “quasi” constant term (“quasi” in the sense that is given by 
the other fields (p and h that are slow-rolling), plus the pure oscillatory term with an amplitude 
that decays with the expansion. So when compared to the evolution of the other fields, (p and 
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h, in a few e-folds the oscillatory term in N will be averaged to zero, and if we require a local 
minimum |^ = 0 then we have effectively, from Eq. (20), 


N{t) ~ 






h\t) 


0(f) , 


( 21 ) 


V2 

which relates the inflationary trajectory of the average value of N to the slowly rolling h, 0 helds. 
Therefore the oscillations of N are not exactly around zero but its centre will move along with 
the inflaton and Higgs held. Therefore in an effective sense the N{t) held will also slowly roll 
along the valley of minima controlled by the Higgs and the inhaton. Ehectively the three helds 
will follow a slow-roll trajectory. 


0 (f) - , (22) 

h{t) ~ -r]hHh{t ), (23) 

N{t)c:^-{r]^ + 2r]h)HN{t), (24) 

with N(t) h(t) 0(f), as in Eq. (17) and rjh = with ? 7 <^ = m?^/{3H‘^) as in Eq. 

(14). 


3 Evolution of the fluctuations during inflation 

During inhation with several light (slow-rolling) scalar helds 0 q,, the comoving curvature pertur¬ 
bation [30] can be written as [35] 

(25) 

in terms of the gauge-invariant scalar held amplitude perturbations, the Sasaki-Mukhanov vari¬ 
ables Qi [31], 

Qa = h4>a + ^0- (26) 

And for the entropy perturbations [36] we have^: 

~ -3i^(0 - % , (27) 

<Pa 0/3 

These expressions may be compared to the (more intuitive) expressions given earlier in Eqs. 
(1), (5), whereTT-a = HQa/cf). More general expression can be found in Appendix B, Eqs. (65) 


^This expression for SafS in terms of the scalar field perturbations holds as far as we can neglect the coupling 
between the scalar, and the ratio Qaj<t>a — Constant, both of which are a good approximation during inflation. 



and (61). The initial curvature perturbation, Eq. (25), will be dominated by the held with the 
largest velocity, whilst the held with the smallest velocity dominates the entropy perturbations 
in Eq. (27) [12,32-35]. 

At hrst sight one might think that the N held, with the smallest velocity N, according 
to Eq. (24), would give the largest isocurvature perturbation. However, since its evolution is 
controlled by the Higgs held, the ratio Qn/N turns out to be of the same order of magnitude as 
Qh/h, as discussed in Appendix A. Thus the Higgs helds contribute strongly to the isocurvature 
perturbations, and in addition dominates the evolution of the N held during inhation according 
to Eq. (21). 

During inhation the curvature perturbation is dominated by the held with the highest veloc¬ 
ity, which according to Eq. (22) is the inhaton. For a given quantity A, the spectrum is dehned 
by [1], 

P 

P^{k) ^—{\A\^), (28) 

where k is the comoving wavenumber. Using Eq. (25) we hnd 


pl/2 

-'7^ 


27r 


H, 


(29) 


where the subscript denotes the time of horizon exit, H^a = k. This is by far smaller than 
the required COBE value ~ 5 x 10“^ [8], unless we take Km^ ~ 10“^® GeV, i.e. an inhaton 
mass of the order of a few eV. Such a tiny value^ is far smaller than the MeV value we need 
to satisfy the slow roll conditions, as discussed previously. Therefore it is natural to suppose 
that the curvature perturbations are initially too small to satisfy the COBE condition, but 
are generated by the conversion of the isocurvature perturbations from the Higgs held, during 
reheating, as discussed in the next section. 

After horizon exit, the entropy and curvature perturbations evolve independently until the 
end of inhation. Following Ref. [12], by the time inhation ends the amplitude of the huctuations 
can be given in terms of their values at Hubble crossing (aif* = k), 


n\ 


H, 


H, 


S4>h\i ~ -y—Q/i* — 
3 h* 


8 Q, 


3 rjhK ’ 


(30) 

(31) 


^The tiny mass could be due purely to radiative correction, 
to zero [26]. 


P{K(j)c)'^, if the soft-breaking mass is set 
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with 

^ vw ■ 

Clearly, given that h 0, we have IZ S^h |i, so that the isocurvature perturbation from the 
Higgs is nice and large, so that when it gets converted into the curvature perturbation during 
reheating it can account for the COBE observation, as we now discuss. 

4 Evolution of the fluctuations during reheating 

The slow-rolling regime ends when the inflaton held reaches the critical value 0c- At this point, 
the effective mass of the N held changes sign, and it starts moving towards the global minimum of 
the potential at Nq ~ O(0c). The ehective masses of 0 and h also start increasing, and they also 
move quickly towards their global minimum values. Around the global minimum, the ehective 
masses rhj for all the helds are of the same order of magnitude, rh^ c:=l ffiN — fhh ~ ~ 

0{lTeV). Therefore, they all start oscillating with similar frequencies and amplitudes of the 
order of 0c. In this way, the vacuum energy dominating during inhation is equally redistributed 
among the three helds, ~ Pat ~ p/i, with the energy density of the oscillating helds behaving 
as matter. 

This will be the situation until the helds decay and transfer their energy density to radiation. 
Due to the smallness of the singlet couplings, they are very long-lived, with a decay rate of the 
order Tat^^ ~ ~ 0(10“^^GeV), where k ~ A ~ 10“^° whose smallness is explained in 

[25]. Reheating lasts until the singlets completely decay, roughly around the time ~ T,^. On 
the other hand, the Higgs held decays much faster through its gauge interactions, with T^ ~ 
awfdh ~ 0(10 GeH). Because of this, the decay of the Higgses can be considered as practically 
instantaneous when compared to that of the singlets. However, given that Vh/ffih ~ aw ~ 0.01, 
they will still have time to oscillate several times before decaying. Notice also that Ta 3> at 
the beginning of the oscillating phase. From that point of view, the decay of the Higgs can also be 
considered “instantaneous”: the Higgs decay before the energy densities in the oscillating helds 
have been redshifted, so that through the decay a fraction of the vacuum energy dominating 
during inhation is transferred into radiation. Given that during the decay we can neglect the 
ehect of the expansion, we will also assume that the huctuations in the Higgs held at the end of 
inhation are converted into those of the radiation huid, with S^r = S^h\i as the initial condition 
during reheating. 

Therefore, during the reheating era, we are left with a two component huid, made of a 
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mixture of matter (0 and N oscillating) and radiation pr (Higgs decay products), but such 
that initially ^^(0) ~ Pr(0). The evolution equations for the system can be written as 


p^ + SHp^ + T^p^ — 0 , (33) 

Pr + 4:Hpr -T^p^ = 0. (34) 

It will take a while for the initial large fraction of radiation to feel the effect of the inflaton decay 
products, and at the beginning it is redshifted as usual like pr oc with p^ oc a~^ and a being 
the scale factor. It is only later when the contribution of the decay products become comparable 
to the already present radiation that the system is fully coupled, and we have pr oc The 

decay is completed by the time ~ T^, and we are left again with radiation which redshifts in 
the usual way pr oc a~'^. 

For comoving scales k well outside the Hubble radius, k aH, the evolution equations for 
the fluctuations during reheating can be approximated by^ by 

din a {Sp^ + Apr^ 

dS(j)r \‘2pr „ T,^ p^ / 3P(^ P''" \ Q 


v 2 ^ 4 ’'!' ’ 


dS(pr 
din a 


dlna Sp^ + ipr H Pr \4pr + 

Thus, at the beginning of reheating when pr p^, the large isocurvature perturbations coming 
from the Higgs helds in Eq. (31) which dominates the entropy perturbations, acts as a strong 
source for the curvature perturbations IZ. Thus the curvature perturbations IZ quickly grow 
to become of the same order of magnitude as the initial entropy perturbation, and afterwards 
remains practically constant. 

The transfer from entropy perturbations to curvature perturbations is illustrated in Fig. (1), 
based on a numerical integration of the exact evolution equations for the perturbations given in 
Appendix C, Eqs. (78-80), together with Eqs. (33-34), with the boundary conditions in Eqs. 
(30) and (31). From the numerical integration, we hnd asymptotically, 

TZ ~ I . (37) 

5 

Hence, using Eqs. (31), (47), we end with the spectrnm of adiabatic perturbations 

^ Ipi/s--^-^ , (38) 

_ ^ 5 15pfe27rh*’ ^ ^ 

^The exact equations for a two-component fluid obtained from Ref. [36] are given in Appendix B for com¬ 
pleteness. They also include the equations for the relative velocity perturbations Vap. In Appendix C we give 
the equations and approximate solutions for the scenario we are considering. In this case, in order to understand 
the qualitative behaviour of the numerical solution is enough to consider the evolution equations for TZ and Safi, 
without including Fi/j, Eqs. (81) and (82) in Appendix C. 
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Log,„(a/a„) 

Figure 1: In the upper panel we show the evolution of the curvature TZ (solid line) and entropy 5'^^ (dashed 
line) perturbations during reheating. Both values are normalised to the initial value of We have taken 

= lO^GeV, (j)/h = 10^° and rj^ = rju ■ In the inset in the upper panel we show a 
blow-up of the initial region on a linear vertical scale. In the lower panel we show the evolution of the energy 
densities and pr, normalised to their initial values. 

which for r]h ~ 0.03, ~ lO^^GeV^ and 27r/t* ~ ITeV^ gives the right order of magnitude as 

measured by COBE. 

The entropy perturbation also remains practically constant until the terms due to the decay 
of the singlets N and 0 become relevant, and the radiation starts behaving as pr oc the 

same as the ratio P(i)/H. From Eq. (36), the evolution for S^r is now approximately given by: 

Sff,r d In a P0 H Pr 

The main contribution in the above equation comes from the second term, which remains con¬ 
stant in this regime with V^p^/{Hp^) — 2.5, and then S^r oc Therefore, when the system 

of matter and radiation becomes effectively coupled, the entropy perturbation actually decreases 
due to the relative increase of radiation coming from matter, and this effect is also seen in Fig. 
(1). Obviously, in this simplest model of reheating once the inflaton completely decays and we 
are left only with radiation, the relative entropy perturbation vanishes in any case. We have 
just integrated the equations only until the time F^ ~ id, without continuing the evolution 
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afterwards when all the energy density in matter becomes radiation. However, it is clear that 
the value of the curvature perturbation, which becomes practically constant near the beginning 
of reheating, will not subsequently be affected by the details near the end of reheating and 
subsequent evolution. 

5 Discussion 

It is worth pointing out a few subtleties regarding details of the transition from inflation to re¬ 
heating. A detailed description of such transition may be relevant for a more accurate estimation 
of the hnal value of TZ. Here we raise some of the concerns related to this point, and point out 
to what extent our conclusions can be affected. During inflation the background evolution equa¬ 
tions for the helds are approximately given by the slow-roll equations, Eqs. (22)-(24); whereas 
reheating starts with the helds already oscillating and their energy densities averaged like matter 
until they decay into radiation. During the transition from one to another, the kinetic energies 
of the helds start increasing and they cannot be neglected any longer. Indeed it is going to be 
the N held which hrst starts moving away from the false vacuum once 0 reaches the critical 
value, with its held-dependent squared mass becoming negative, and its vev and velocity growing 
accordingly. Its quantum huctuations Qn will also feel the same instability, and increase until 
the ehective mass become positive again and the oscillations begin. Given Eq. (25), this may 
mean that during the transition period it is the N held, with both its kinetic energy and the 
amplitude of its huctuations increasing relatively to the others, which will end dominating the 
curvature perturbation. Still, our conclusion that the Higgs is responsible for the hnal value of 
the adiabatic perturbations holds, although indirectly. It is through the coupling of the Higgses 
to the N held that the quantum huctuations of the latter are not suppressed during inhation, 
as it should correspond to a massive held (see Appendix A). We also note that — Qh/h, 

so that even if the N kinetic energy dominates before the oscillations begin, we will end in any 
case with IZ ~ HQj^/N ~ HQ^/h, and our order of magnitude estimation holds. However, 
this assumes that the ratio Q^/N remains constant once the held N start moving toward the 
global minimum. This does not seem unreasonable given that both the background held and its 
quantum huctuations will feel the same instability in the ehective mass. However, it may happen 
that as the helds N and h move toward the global minimum with increasing speed N and h, 
their quantum huctuations still behave as those of a light held without changing appreciably. 
This will tend to decrease the hnal ratios Qn/N and Qh/h, and therefore the value of IZ. 
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On the other hand, we have started evolving the perturbations during reheating after the 
Higgs held decay into radiation. As mentioned in the previous section, before the Higgs decay 
the helds still have time to oscillate several times, such that the energy density is redistributed 
among the three oscillating scalar helds. During the initial stage of the oscillations, as far as the 
three species behave all on average like matter, the curvature perturbation TZ remains practically 
constant on super-Hubble scales and it is not ahected by entropy perturbations (see Eqs. (63) 
and (75) in Appendix B). In other words, the change in the curvature perturbation is given 
by the large-scale non-adiabatic (entropic) pressure perturbation, which is negligible as far as 
the scalar helds behave with similar ehective equation of state. However, this argument does 
not take into account the possible parametric amplihcation of the scalar quantum huctuations 
in a background of oscillating helds, a process that in the context of inhation generically is 
known as “preheating” [37,38]. Through parametric resonance, induced by the time dependent 
ehective mass term in the evolution equations, held mode amplitudes can grow exponentially 
with time within certain resonance bands in k-space. Thus, preheating become a more ehective 
mechanism of transferring energy from the oscillating background helds to quantum huctuations 
than the standard perturbative decay of the helds. The question then is whether it also provides 
a diherent and efficient source for non-adiabaticity, with super-Hubble {k aH) huctuations 
exponentially amplihed during preheating [39,40]. This would translate into an exponential 
increase in the curvature perturbation which would be somehow difficult to keep at the level of the 
COBE observational value. As usual, the answer is model dependent. For example, it has been 
pointed out that for a massive (mass larger than the Hubble parameter) held during inhation 
the mechanism is inefficient, because previous to preheating huctuations are hrst exponentially 
damped during inhation [41]. However, this restriction does not apply to light helds during 
inhation [42], as in the model presented in this paper. Moreover, in hybrid inhation the ehect 
can be stronger due to the presence of negative ehective squared masses at the end of inhation 
[43,44], which are per se a source of instabilities in the evolution equations for the huctuations. 
To which extend the combined ehect on super-Hubble scales is relevant is then a question of 
when the resonance ends [45]. Sooner or later, this occurs due to the backreaction ehect of the 
small-scale (sub-Hubble) modes which also grow resonantly during preheating; in some models 
it may happen before the curvature perturbation exceeds the COBE value, depending on the 
value of the inhaton coupling constant. Again, the issue is highly model dependent. Preheating 
in the present model has been studied previously in Ref. [43], but without neither taking into 
account the Higgs huctuations nor metric perturbations. It was shown that the combination of 


14 



parametric resonance plus tachyonic instability gives rise to the exponential grow of the inflaton 
and N held huctuations, and the resonance was stronger on the small-scale perturbations around 
k ~ 0.3m^. At the same time and because of that, backreaction effects become important very 
quickly after only three oscillations of the helds. Given that the tachyonic instability will be 
still present mainly in the effective N mass, we do not expect this conclusion to change once the 
Higgs huctuations are also taken into account. Even they may speed up the end of the resonance, 
keeping under control the growth of the curvature (and entropy) perturbations. Moreover, an 
increase due to parametric resonance might compensate for a potential decrease during the end 
of the slow-roll regime. Nevertheless, the hnal answer would require the numerical integration of 
the evolution equations, for both the background helds and quantum huctuations, from inhation 
to reheating, which is beyond the scope of the present publication. 

6 Spectral index, gaussianity and gravity waves 

In this section we discuss the main predictions of the mechanism discussed here as compared 
to both the standard curvaton mechanism, based on the late-decaying scalar, or on standard 
hybrid inhation. 

In the original hybrid inhation model based on the NMSSM [20], where the inhaton was 
responsible for curvature perturbations, we predicted a hat spectrum of curvature perturbations 
with a spectral index indistinguishable from n = 1. In the present NMSSM hybrid inhation 
model, in which the Higgs has large isocurvature perturbations which get transferred to curvature 
perturbations during reheating, we would expect a spectral index controlled by the Higgs helds 
instead of the inhaton and given by^ 

n - 1 = 2rih - ^ ‘2r]h ~ 2m\/ (40) 

The prediction for the spectral index in Eq. (40) involves the soft mass parameter ruh which 
depends on the details of the supersymmetry breaking mechanism. Eor example in the extra- 
dimensional model introduced in the companion paper [25], Higgs soft masses are not expected 
to be much suppressed with respect to the Hubble rate of expansion [25]. Using n = 0.93 ± 0.13 
[2] we have an upper bound on the spectral index n < 1.06 which leads to a constraint on 
the soft Higgs mass of rjh < 0.03. The inhaton mass on the other hand is only constrained by 
the slow-roll condition tj^ < 1, since the inhaton does not contribute signihcantly to curvature 

^This prediction for the spectral index is valid when there is no relic isocurvature perturbation which can 
significantly affect the CMB power spectrum. 
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perturbations. We would therefore expect typical deviations such as | n — 1 |~ 0.1. If n = 1 
is measured very accurately then this would be evidence for the original inflaton generated 
curvature perturbations. 

In the original single field inflaton model [20] we predicted a Gaussian CMB spectrum. In 
the present case we would expect small non-Gaussian effects, as discussed in Ref. [16], which 
contribute quadratically to the Higgs energy density. 


5ph Jh , {6hf 

— ~ z -- 

Ph h 


(41) 


Given that the amplitude of the held perturbations at horizon crossing are of the order of the 
Hubble parameter, (5/i ~ iJ* ~ 10 MeV, but the typical value for the vev of the Higgs is h ~ 1 
TeV, the linear term clearly dominates in Eq. (41) given a Gaussian spectrum. The small 
non-Gaussian effects can be parametrised by [46] 


IInl 



SPh/Ph 

n 



(42) 


where we have used the approximations in Eqs. (31), (37) and (41). Within the approximations 
involved in the above equations, and using ph < 0.03 we then get | Jnl \< 0.3, which is below 
the expected upper bound by the PLANGK satellite | Jnl \< 5 [46]. Therefore we do not expect 
signihcant (observable) non-Gaussian effects in this model. 

Goncerning gravity waves, as is typical of hybrid inhation models with a low inhation scale, 
the previous model [20] predicted negligible effects in the GMB spectrum from gravity waves. In 
the present approach we similarly expect small effects of gravity waves in the GMB spectrum. 
Gravitational waves are generated with a spectrum [1] 

~ 8^2— , (43) 

nip 

and the tensor-scalar ratio, from Eqs (37) and above, is then given by 

^ ~ {UOTTfeh oc f—) , (44) 

where the value of the Higgs held at the time of horizon exit is typically h* ~ ITeH, which is 
negligible compared to the reduced Planck mass mp. Therefore we do not expect any observable 
effects of gravity waves in the GMB spectrum. 

It has been pointed out that in general for the curvaton models based on a late-decaying scalar 
the entropy perturbation originated during inhation might survive the reheating era as a relic 
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isocurvature perturbation between radiation and some other present component, say cold dark 
matter (neutralinos) or neutrinos [15,16]. In the framework described in this paper the Higgs 
helds can decay into neutralinos and in addition give rise to leptogenesis during the reheating 
process [47]. Therefore the Higgs isocurvature perturbations might be expected to give rise to 
relic isocurvature perturbations both in neutralino dark matter and in baryons, as compared to 
photons. However the Higgs excitations decay near the beginning of the reheating era, and any 
neutralinos produced at that time are relativistic and share the curvature perturbation with that 
of the photons. Therefore once the neutralinos decouple from the radiation fluid, their curvature 
perturbation will remain constant and equal to that of the radiation at the time. On the other 
hand throughout reheating the radiation will be coupled to the matter energy density through 
the decay of the inflaton field. For such a coupled system, individual curvature perturbations 
are not conserved, unlike the late-decaying scalar scenario in [10,11,14-17], and so in our case we 
expect the isocurvature perturbation of the photons to change in a complicated way during the 
reheating process. Thus in our case the hnal neutralino isocurvature perturbation will depend on 
the details of the reheating process, in particular when the neutralino decouples and becomes non 
relativistic during the reheating process. It would be interesting to explore this quantitatively 
in a future work. 


7 Summary 


To summarise, we have proposed and discussed a new implementation of the curvaton scenario 
[10,11,14-17] in which quantum fluctuations of light scalar helds other than the inhaton can be 
responsible for the generation of the primordial anisotropies in the Universe. Unlike the previous 
scenarios, we do not assume a late-decaying scalar which is decoupled from the inhaton held. 
Instead we have proposed a new mechanism based on supersymmetric hybrid inhation in which 
the isocurvature perturbations of a hybrid held coupled to the inhaton could be transferred to 
curvature perturbations during the initial stages of reheating. We have further suggested that 
good candidates for such helds are the Higgs doublets of the supersymmetric standard model, 
which may have a hat direction during inhation provided that the fi term is generated after 
inhation. We have discussed in detail a specihc model which implements this mechanism, based 
on a supersymmetric hybrid inhation model which is a variant of the NMSSM, based on the 


superpotential in Eq.(6). We have shown that for this particular model our mechanism leads to 


pl/2 

-' 7 ^ 


rsj 


H, 

27rh* ’ 


(45) 
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which for a typical value of the Higgs vev 27rh* ~ 1T eV gives the correct order of magnitude as 
measured by COBE 10 This is an important success of both the particular model and 

of the proposed mechanism in general. The main prediction is that the spectral index is expected 
to deviate significantly from unity | n — 1 |~0.1. Ifn = 1 is measured very accurately then this 
would be evidence for the original inflaton generated curvature perturbations. No observable 
signals of non-Gaussianity, or gravity waves are expected in the CMB spectrum. However there 
may be observable relic isocurvature perturbations between radiation and some other present 
component, say cold dark matter (neutralinos) or neutrinos which depend on the details of 
reheating and differ signihcantly from the expectations of a late-decaying scalar [15,16]. 

To conclude, we have proposed and discussed the attractive possibility of having the Higgs 
held of the supersymmetric standard model as being responsible for the large-scale structure 
of the Universe, providing a further strong link between cosmology and particle physics of 
the kind recently emphasised in [48]. Eventually, if such a theory as that presented here is 
realised in nature, it should be possible, from both laboratory and cosmological measurements, 
to demonstrate the links between particle physics and cosmology inherent in such a model [48]. 
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Appendix A 


The equation of motion for the gauge invariant quantum huctuation Qq, with comoving wavenum¬ 
ber k, can be written as [49] 


Qa + SHQa H-jQo + 


V^P- 


a^m% V H 


Q/3 — 0, 


(46) 


where = d'^V/d<pad(j)i 3 . For the inflaton and Higgs helds during inflation, we have Vhh, 
V^h ^ and therefore the fluctuations and Qh will be frozen to a constant value once 
outside the horizon, k < aH, given approximately by the value at horizon crossing, 

" TIP ’ 
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for a = (j), h. On the other hand, neglecting for simplicity the sub-dominant terms coming from 

metric contributions, the evolution equation for Qtv reads 

^2 

Qn + 3HQ]\f -|- (— -|- Vnm)Qn + VN<j>Q(j> + VNhQh — 0 . (48) 

CL 


Like in the case of the evolution of the background held N, now the large mass term Vnn ~ 
0{k,‘^(PI) 3> gives rise to oscillations with an amplitude decaying as a~^, but displaced from 
zero due to the and Qh terms. Therefore, after horizon exit, and averaging over the fast 
oscillations in a Hubble time, the N held huctuation will also tend to a constant value given by 

Qk-^Q4. + ^Qi.- ( 49 ) 

^NN ^NN 

Using Vnn — '^n 4 > — — Ak/i^, V]\fh — 2A/t0ch, we have 


Qn 


Therefore, the huctuation Qn is 
background held N is suppressed 


^ ( 2 - 


4>c 


\h ^, ^ \h ^ 

2)^0 _ 


n4>c 


Xh 

K(j)c 


Qh 


(50) 


suppressed with respect to Qh by the same factor than the 
with respect to h, Eq. (21), such that. 


and then 


TV ~ - 

Xh ■ 

n4>c 

(51) 

Qn 

Qh 

(52) 

_ 

N 

h ■ 


Appendix B 


In this Appendix we summarise our conventions and notation for the perturbations, and give the 
evolution equations for curvature and entropy perturbations in a multi-component huid. These 
can be found in Refs. [36,12]. Using linear perturbation theory, the equations are given in terms 
of gauge-invariant quantities Ac = 5pc/p, V and p describing the amplitude perturbation of 
the total density in the comoving frame, velocity and entropy respectively. In addition we have 
the relative gauge-invariant variables between any two components a and (3: entropy and 
relative velocity Vap- 

Gauge-invariant quantities are dehned from the original perturbations in the stress-energy 
tensor and the metric. Linear scalar perturbations of the metric are given by the line element: 

ds'^ = —(1 -I- 2A)dt^ + 2aBidx^dt + a^[(l — 2hp)5ij + 2Eij]dx^dx^ , (53) 
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where ifj is the gauge-dependent metric perturbation. Combining B and E, one can dehne 
the gauge-invariant variable 

$ = - Ha{B - aE ), (54) 

which is the curvature perturbation in the longitudinal (or zero-shear) gauge {aEi — 5; = 0). 

Each fluid component is described by a perfect-fluid stress-energy tensor, with background 
energy and pressure Pa- Although the total stress-energy tensor is conserved, the stress- 
energy tensor of each component may not be conserved individnally. Therefore, the nnperturbed 
continuity equation for a given component a is given in general by: 

Pa + SH(pa Pa) = Ca ; (55) 


where a dot denotes derivate with respect to time, and Ca are sonrce terms subject to the 
total energy-momentum conservation constraint J2a Cq, = 0. For example, in Eqs. (33) and 
(34) we have = —Cr = —T^p^. For later use, we define the equation of state for each 
component, Wa = Pa/Pa, and the sound velocity = Pa/Pa- To simplify notation, we also 
define ha = Pa + Pa- 

Pertnrbations in the stress-energy tensor are given by 5pa, 5Pa and the momentnm pertur¬ 
bation Sqa (neglecting the anisotropic stress). Gauge-invariant variables Spca, To and rja are 
defined as: 


snch that 


hpca 

Va 

PaPa 


— hpa + P, 

k ( 


5q 


V + P’ 

he 


n \PaP Pa 

hPa ^sa^Pci , 


— aB + a E 


hpc 

ip + P)V 

Sap 

Vap 


'y h Pea , 

a 

T Pa)Va , 

a. 

h Pea h PeP 


Pa + Pq Pp P Pp’’ 


= Va-Vn- 


(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


The total entropy pertnrbation (or non-adiabatic pressnre perturbation) Pp = 6P — c^Sp can 
be written in terms of the individnal entropy perturbations pa for each component and Sap as: 


w 


1 + w 


V = J2 


he 


WaPo 


h 1 + Wo 


1 hahp ^^2 

2^2^ h2 

a p 


(cL ~ <^lp)Sap + 


A. 


1 +W 


E' 


a 

‘3Hh 


(63) 
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Individual entropy perturbations rja vanish for adiabatic perturbations, but cannot be neglected 
for example in the case of a scalar held, in which case they are given by 


PaVci (1 ^sa)i.^Pca j Pai^a ^)) 


The comoving curvature perturbation is dehned by: 

where in the last equality we have used the relation between $ and the comoving total density 
perturbation through the energy constraint given by 

3 faH\‘^ . 


2\k J ^ 

On the other hand, the curvature perturbation on constant-density hypersurfaces is given by 

C = - H— , (ei 

P 

which is related to the comoving curvature perturbation IZ by 


-C = n- 


9{p + P) \aH 


Therefore, for super-Hubble scales, k aH, ( and IZ are approximately the same. 

We note that we have followed Refs. [12,36] in dehning the entropy perturbation Sap. An 
alternative dehnition related to the individual curvature perturbations Ca, used in Refs. [10,16], 
is given by 

= -3hf = 3(Ca - 0) • (69) 

\ Pa PP / 

In terms of Sap, the total entropy perturbation is then given by: 


1 + w 




hoL '^a^a , 1 ho^hg .9 2 \ n 


Nevertheless, for a system of uncoupled fluids, Eqs. (61) and (69) become identical. Sap and 
Sap differ by terms proportional to the source terms Ca when there is energy transfer between 
components. For example, for a 2-component system, their relation can be written as 

(1 - Qa){l - qp)SaP = Sap + {qa “ qp)^^ > ( 71 ) 

where 1 — qa = —pa/{3Hha)- 
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The evolution equations for the total density and velocity fluctuations in a flat Universe are 
given by: 


= 3tcAc - (1 + w)-^U, 
an. 


= -V-l( 


v = -V 


3 /all 


2\k 


Ac + 


k 


aH \ 1 + tc 


Ac + 


w 


1 + tc 


V 


(72) 

(73) 


where prime means derivative with respect to In a. Thus, taking the derivative in Eq. (65) and 
using Eqs. (72-73), the evolution equation for the comoving curvatnre pertnrbation reads 

2 o / „ rr \ 2 


/ n3 faH\ . 3 faH\ ^ 

K = (l + 3«.)- — A,-- — + + 


2 \ k J ° 2 \ k J ° 2' ' k 

ct . w ‘2 I k w 

1 + w ^ 1 + w^' 


k 


^^Ac + 

1+w 1+w 


(74) 


(75) 


Finally, for a two-component flnid, snch that Ca = —C/ 3 , the equations for the relative 
fluctuations Saf 3 and Va /3 are given by: 


c> _ 
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Vap — 3 


WaVo 


wprjp 
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1+Wa 1+Wp^ 
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(1 + ^sa)j^ - (1 + ^sp'^J^ ) > 


(77) 


where Eca = Eop are the gange invariant pertnrbations of the energy transfer term in the 
comoving frame. For the simplest case Ca = TaPa then Eca = hpca/Pa- 


Appendix C 

For the scenario considered in this paper, we can take the component “a” to be the oscillating 
fields behaving like matter {wa = Csa=0),“/5” then refers to the radiation {wp = Csp=l/‘i) 
initially coming from the Higgs decay products, and Ca = Cp = —Tppp. Then, the evolution 
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Eqs. (75), (76) and (77) for the curvature 71, and V^r = {,Ha/k)V^ reduce to: 
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“ 2 '^^^^ + 3p^ + 4pc^^^ + 31 + u; + 3p^ + 4pc^^^ 

, r, ( 3pl-3pl \ 

H \pr{3p^ + Apr)] 

For super-Hubble scales such that k aH, we can further simplify the system by neglecting the 
contribution of the terms oc {k/aH^, and that of Ac = 2{k/aH)^^/3. In order to understand 
the qualitative behaviour of the solutions to the above system of equations, we can focus on 
the evolution of the curvature 7Z and S^r, without taking into account the terms due to V^r- 
From the evolution equations we can see that S^r and V^r will behave approximately the same, 
so both of them will affect the curvature 77 in a similar way; then the main effect is obtained 
by keeping only those term due to S'^c in the evolution equations. We have checked that the 
numerical solution of Fqs. (78-80) is well approximated by that of: 

77' ~ -^Mj-—Ssr (81) 

(3p^ + 3pc)2 ^ ^ 

o' _. 12pc „ p^ / 3p0 + Pr ^ r, 


77' ~ 


Q' _9 , 

■ 3p^ + 4pc^^^ 


P<j> ( ^Ptp "h Pr \ (~I 

H Pr V3p0 + 4pc/ ’ 


In order to solve analytically the equations, we can distinguish two regimes depending on the 
evolution of the background energy densities p^ and pc, as it can be seen in Fig. (1): (a) from the 
beginning of the reheating period at oq, with pc(ao) — p^^ao), up to say ai, the decay products 
coming from the singlets have no effect on the radiation, and both components, matter and 
radiation, behave as if they were decoupled, with pr oc and p^ oc a“^; (b) from ai to the 
end of reheating, the singlets decay products start contributing to the background radiation, 
with Pr oc and p^ oc a~^. Therefore, for Oq < a < Oi, Fqs. (81) and (82) can be solved 

neglecting the terms proportional to V^/H, and we obtain for 77 and S^p- 


77 ~ 77(ao) + — V(«o) ( 5 
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where 7^(ao) 5 ' 0 r(ao) are the initial values of the perturbations taken as those at the end of 

inflation. It can be seen that very quickly TZ will tend to a constant value with TZ ~ S' 0 r(ao)/lO. 
Similarly, the entropy perturbations go to ~ 0.08 S'<^r(0). That is, entropy perturbations act 
as a source for the curvature perturbations only at the beginning, as far as the ratio Pr/P 4 > is 
not suppressed. Afterwards, pr p^p, and 7Z remains constant. 

Entropy perturbations will not further affect the evolution of the curvature even when the 
background equations become coupled at Oi, and the effects of T^ cannot be neglected. At 
this point we can still assume that the energy density is dominated by that of the singlets, 
behaving like matter, with H oc Therefore, from oi onwards, the combination V^ppp,/[Hp^.) 

remains constant until the singlets completely decay. Numerically we hnd Vp,pp,/{Hpr) — 2.5, 
independently of the value of Taking this into account, we can now solve Eq. (82) including 
the term: 

S^{a > ai) ~ 0.07 S^r{ai) . (85) 


That is, the relative perturbation between matter (singlets) and radiation decreases. With this 
result, we can see that the source term in the evolution equation for the curvature perturbation, 

TZ' ~ oc ] , ( 86 ) 

3p0 Vai/ 

also decreases in time and it is no longer effective, even if by the end of the reheating period 
Pr ~ P 4 , again. 

The evolution of the perturbations plotted in Fig. (1) is the result of the numerical integration 
of the system of Eqs. (78-79), with initial conditions at the end of inflation [12]: 
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( 88 ) 
(89) 


Given that for the amplitude of the held perturbations at the time of horizon crossing we have 
— Qh* — 77*/v^2^, it is more convenient to normalise the above perturbations with respect 
to the initial value of the entropy perturbation Sp,r\i-, such that: 

n 
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